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1. Introduction and motivation 

Although topologically nontrivial configurations of the gauge fields allow for a CP-violating 
term in the Lagrangian of QCD |2p, experiments indicate that its numerical coefficient, known 
as 6, is vanishingly small, 6 < 10~ 10 ^]. Nevertheless, in spite of the fact that spontaneous 
breaking of P and CP were proved to be forbidden in the true vacuum of QCD for 6 = there 
is still room for metastable CP- violating states at finite temperature This provides the exciting 
possibility of probing the topological structure of QCD, and leads to the idea that P- and CP- 
odd domains could be produced within the quark-gluon plasma presumably formed in high-energy 
heavy ion collisions []7J]. 

The experimental signatures that were proposed to identify the presence of CP violation in 
certain domains within the plasma rely on charge separation of hadronic matter [^]. This asymme- 
try can be enhanced to the point of becoming detectable by current experiments at RHIC by the 
presence of a very strong magnetic background, as in the case of noncentral collisions, generating 
a rich phenomenon denominated chiral magnetic effect [g]. 

The quark-gluon plasma possibly possesses regions with nonzero winding number, Q w ^ 0, 
where sphaleron transitions are induced. In the particular case of noncentral collisions a strong 
magnetic field is generated in the orbital angular momentum direction, perpendicular to the reaction 
plane. In this process, the strong magnetic field B restricts the quarks (all in the lowest Landau 
level, aligned with B) to move along its direction. For a topologically nontrivial domain with e.g. 
<2m = — 1, left-handed quarks are converted into right-handed ones, inducing an inversion of the 
direction of momentum and, consequently, a net current and a charge difference are created along 
the direction of the magnetic field. The system is, then, P-odd. 

The magnetic fields involved, although short-lived for very high energies [§], will certainly 
live longer for lower values of sfs, and can attain enormous intensities, above those considered for 
magnetars Jl0| ] and comparable only to the ones believed to be present during the early stages of 
the universe [|TI|]. The presence of such extreme magnetic fields coupled with strong CP violation 
raises several theoretical questions. The first one is how the QCD phase diagram is altered by the 
presence of a nonzero uniform magnetic field that plays the role of another "control parameter". 
Next one can ask where are the possible metastable CP-odd states and how "stable" they are, i.e. 
how long their lifetimes are. These questions are directly related to possible modifications in the 
nature of the phase transitions of strongly interacting matter and to the relevant time scales for phase 
conversion. As we have mentioned above, there is already an on-going experimental investigation 



that seems to show indications of the chiral magnetic effect Q12|]. Furthermore, pioneer lattice 
studies on the chiral magnetization of the non-Abelian vacuum also seem to be able to capture 
these effects [||]. 

In what follows, we discuss the effects of a strong and constant magnetic background and of 
CP violation on the chiral transition at finite temperature and vanishing chemical potential. For this 
purpose, we adopt the linear sigma model coupled with two flavors of quarks as our effective field 



theory Q14], following the notation and conventions of Ref. [ 150. We show that for high enough 



magnetic fields the chiral transition is no longer a crossover. Instead, it is turned into a first-order 



transition [16]. 



To include the effects of the presence of the axial anomaly and CP violation, we add a term 
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that mimics the presence of nontrivial gauge field configurations, the 't Hooft determinant [jT 
The rich vacuum structure brought about by a nonzero 6 term in the action has an influence on 
the chiral transition, and generates a more complex picture in the analysis of the phase diagram 



of strong interactions Q18|]. Working in a mean field approximation in an extended CP-odd linear 



sigma model, our results can be cast in terms of condensates of the fields Q19|]. In this framework, 
three different phases occur, making the topography of extrema very rich, and allowing for the 
existence of metastable minima in certain situations. 

2. Effective theory for the chiral transition in a strong magnetic background 

Modifications in the vacuum of CP-symmetric QCD by the presence of a magnetic field have 



been investigated previously within different frameworks, mainly using effective models [2C, 21 



2|, H ||, H, ||, 0], especially the NJL model [gf], and chiral perturbation theory |(| ||], 



but also resorting to the quark model [ |32| ] and certain limits of QCD Q33|]. Most treatments have 
been concerned with vacuum modifications by the magnetic field, though medium effects were 
considered in a few cases, as e.g. in the study of the stability of quark droplets under the influence 
of a magnetic field at finite density and zero temperature, with nontrivial effects on the order of the 



chiral transition []34[]. More recently, magnetic effects on the dynamical quark mass [35] and on 



the thermal quark-hadron transition J36Q, as well as magnetized chiral condensates in a holographic 



description of chiral symmetry breaking [37], were also considered. Recent applications of quark 



matter under strong magnetic fields to the physics of magnetars using the NJL model can be found 



in g. 

To investigate the effects of a strong magnetic background on the nature and dynamics of the 
chiral phase transition at finite temperature, T, and vanishing chemical potential, we adopt the 



linear sigma model coupled to two favors of quarks. Following the notation of Ref. Q15|], we have 
the lagrangian 

S£ = Wf[ifd^-g(<y + iY 5 r-7i)] Yf + ^idfiOd^o + d^^-Vio,^) , (2.1) 



where V (a, k) = ^(o 2 + K 2 —v 2 ) 2 — ho is the self-interaction potential for the mesons, exhibiting 
both spontaneous and explicit breaking of chiral symmetry. The Nf = 2 massive fermion fields 
\j/f represent the up and down constituent-quark fields y = (u,d). The scalar field a plays the 
role of an approximate order parameter for the chiral transition, being an exact order parameter for 
massless quarks and pions. The latter are represented by the pseudoscalar field n = (n° ,n + ,n~), 
and it is common to group together these meson fields into an 0(4) chiral field = (a, n). In what 
follows, we implement a simple mean-field treatment with the customary simplifying assumptions, 
where quarks constitute a thermalized fluid that provides a background in which the long wave- 
length modes of the chiral condensate evolve. At T = 0, the model reproduces results from chiral 
perturbation theory for the broken phase vacuum. In this phase, quark degrees of freedom are 
absent (excited only for T > 0). The a field is heavy, M a ~ 600 MeV, and treated classically. 
On the other hand, pions are light, and fluctuations in 7i + and 7i~ couple to the magnetic field, B, 
as will be discussed below, whereas fluctuations in 71° give a B-independent contribution that we 
ignore, for simplicity. For T > 0, quarks are relevant (fast) degrees of freedom and chiral sym- 
metry is approximately restored in the plasma for high enough T. In this case, we incorporate 
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quark thermal fluctuations in the effective potential for a, i.e. we integrate over quarks to one 
loop. Pions become rapidly heavy only after T c and their fluctuations can, in principle, matter since 
they couple to B. The parameters of the lagrangian are chosen such that the effective model re- 
produces correctly the phenomenology of QCD at low energies and in the vacuum, in the absence 
of a magnetic field. Standard integration over the fermionic degrees of freedom to one loop, us- 
ing a classical approximation for the chiral field, gives the effective potential in the a direction 
V e ff = V(<p) + V q (0), where V q represents the thermal contribution from the quarks that acquire an 
effective mass M(a) = g\o\. The net effect of the term V q is correcting the potential for the chiral 
field, approximately restoring chiral symmetry for a critical temperature T c ~ 150 MeV [|l5|]. 
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Figure 1: Evidence for a first-order chiral tran- 
sition in the effective potential for eB = \0m\. 
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Figure 2: Zoom of the barrier for eB = dm % . 



Assuming that the system is now in the presence of a strong magnetic background that is con- 
stant and homogeneous, one can compute the modified effective potential following the procedure 



outlined in Ref. Q16|]. In what follows, we simply sketch some of the main results. For definiteness, 
let us take the direction of the magnetic field as the z-direction, B = Bz. The effective potential 
can be generalized to this case by a simple redefinition of the dispersion relations of the fields in 
the presence of B, using the minimal coupling shift in the gradient and the field equations of mo- 
tion. For this purpose, it is convenient to choose the gauge such that A^ = (A , A) = (0, —By,0,0). 
Decomposing the fields into their Fourier modes, one arrives at eigenvalue equations which have 
the same form as the Schrodinger equation for a harmonic oscillator potential, whose eigenmodes 
correspond to the well-known Landau levels. The latter provide the new dispersion relations 

pl n =p 2 z +m 2 + {2n+l)\q\B , p 2 n = p\ + m 2 + (2n + 1 - o)\q\B , (2.2) 

for scalars and fermions, respectively, n being an integer, q the electric charge, and a the sign of 
the spin. Integrals over four momenta and thermal sum-integrals are modified accordingly, yielding 
sums over the Landau levels. 

In our effective model, the vacuum piece of the potential will be modified by the magnetic 
field through the coupling of the field to charged pions. To one loop, and in the limit of high B, 
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eB » m 2 n , one obtains (ignoring contributions independent of the condensates) [16] 



V\+V v 



2mfcB 



(2.3) 



Thermal corrections are provided by pions and quarks. However, the pion thermal contribution 
as well as part of the quark thermal contribution are exponentially suppressed for high magnetic 
fields, as has been shown in Ref. JT^]. The only part of the quark thermal piece that contributes is 



V 



eBT 

27^ 



J + ~dx In (l+e-V^W 



(2.4) 



where N c = 3 is the number of colors. Therefore, the effective potential is corrected by the contribu- 



tions in (2.3) and (2.4) in the presence of a strong homogeneous magnetic background. Therefore, 
the presence of the magnetic field enhances the value of the chiral condensate and the depth of the 
broken phase minimum of the modified effective potential, a result that is in line with those found 



within different approaches (see, for instance, Refs. pq , |29| , pl|]). 

Fig. 1 displays the effective potential for eB ~ 10m| at different values of the temperature to 
illustrate the phenomenon of chiral symmetry restoration via a first-order transition. For RHIC top 
energies one expects eB ~ 5 — 6m| For lower values of the field, the barrier is smaller. In Fig. 
2, we show a zoom of the effective potential for eB ~ 6m 2 n for a temperature slightly below the 
critical one. This figure highlights the presence of a first-order barrier in the effective potential. For 
a magnetic field of the magnitude that could possibly be found in non-central high-energy heavy 
ion collisions, one moves from a crossover scenario to that of a weak first-order chiral transition, 
with a critical temperature ~ 30% higher [16]. Although the barrier in this case is tiny, the intensity 



of supercooling of the system is expected to be rather large due to the smallness of nucleation rates 
when compared to the expansion time scales for the heavy ion scenario. Therefore, even a small 
barrier can keep part of the system in the false vacuum until the spinodal instability is reached 
and the system is abruptly torn apart. Nevertheless, since the magnetic field falls off very rapidly 
at RHIC top energies, we expect that even for lower values of y/s only the early-time dynamics 
should be affected. 

As caveats, first we note that although non-central heavy ion collisions might show features of 
a first-order transition when contrasted to central collisions, in this comparison finite-size effects 
become important and have to be taken into account [[3^]. Second, since the magnetic field varies 
very rapidly in time, it can induce strong electric fields that could play a relevant role via the 
Schwinger mechanism. 



3. CP-odd linear sigma model 

To describe the chiral phase structure of strong interactions including CP-odd effects, we adopt 
an effective model that reproduces the symmetries of QCD at low energy scales, and has the appro- 
priate degrees of freedom at each scale: the CP-odd linear sigma model coupled with two flavors 
of quarks. The chiral mesonic sector is built including all Lorentz invariant terms allowed by 
symmetry and renormalizability. Following Refs. [|l7|, iO], one can write 



l -Tr(d^d^) + |Tr(^) - ^[Tr(f ^)] 2 - ^Tr[(f ^f] 
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+^[ e ' e det(0) + ^ e det(0 t )]+Tr[/ J ((^ + t )] 



(3.1) 



The potential in the Lagrangian above displays both spontaneous and explicit symmetry breaking, 
the latter being implemented by the term ~ h. The strength of CP violation is contained in the 't 
Hooft determinant term, which encodes the Levi-Civita structure of the axial anomaly and depends 
on the value of the parameter 6. 
Expressing the chiral field as 



(j) = (cr + irj) + (do + Hi) ■ ? 



(3.2) 



where f are the generators of SU(2), the Pauli matrices, the potential takes the following form 
(substituting the parameter h by H = ylM): 



Vy 



--(o +K +T] 2 + 3o) - -cos0 ((7 +K -r\ -afy + c sin0 [at] — K-ao) 

1 X IX 

-Ho + - (Xi + — ■ )(a 2 + T] 2 + 7T 2 + a\f + —{goq + T}% + % x a ) 2 , 



(3.3) 



where the parameters a, c, H, X\ and X2 are fixed by vacuum properties of the mesons [19]. Quarks 
are coupled to the chiral fields in the same fashion as before. 

Following a mean field analysis for the condensates (a) and (tj) and assuming that the re- 
maining condensates vanish, we can compute the effective potential, which is a function of the 



condensates above and of the CP violation coefficient 8, and fix all the free parameters []19[]. For 
6=0 the model, of course, reproduces the results from the usual linear sigma model. Increasing 
the value of 8 from zero to it, the minima of the effective potential rotate from the a direction 
almost to the 77 direction, the rotation being complete only for massless quarks. 
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Figure 3: Contour plot of the effective potential 
for = % and T = 125 MeV. Numerical values 
are in MeV. 



Figure 4: Contour plot of the effective potential 
for = % and T = 128 MeV. Numerical values 
are in MeV. 
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In Figs. 3 and 4 we show contour plots of the effective potential in the case of = n. Increas- 
ing the temperature, the minima move towards the center, indicating chiral symmetry restoration. 
However, there is a clear barrier between the global minimum and the new minimum that becomes 
the true global minimum at high temperature (at tj = 0). In contrast to the case with 6=0, this 
signals a first-order transition, and the possibility of metastable CP-odd states. Moreover, since the 
critical temperatures for the melting of a and tj condensates are different, three different phases are 
allowed in systems with 6 between zero and %: one in which both condensates are present, another 
where the tj condensate vanishes, and a phase where both condensates vanish [19]. This can also 
be illustrated by the behavior of each condensate as a function of the temperature, as shown in Fig. 
5. 
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Following the steps of the previous section, we can incorporate the effects from a strong mag- 



netic background on top of the CP-odd linear sigma model Q19|]- One then finds that the critical 
temperature is raised, as well as the barrier along the tj direction. Effects on the a field are anal- 
ogous to the ones obtained in the CP-conserving linear sigma model. In Fig. 6, we illustrate this 
phenomenon for a magnetic field of eB = 10m| and 6 = %. 

Therefore, for nonzero 6 metastable CP-violating states appear quite naturally in the CP-odd 
linear sigma model. (However, they were not found in an extension of the NJL model [|18|]. For a 
discussion, see Ref. [[fl]].) In fact, larger values of 6 tend to produce a first-order chiral transition 
and might lead to the formation of domains (bubbles) in the plasma that exhibit CP violation JI^]. 
This reinforces the scenario of possible metastable CP-odd states in QCD that are so relevant for 
the chiral magnetic effect [§|. This behavior is enhanced by the presence of a strong magnetic field, 
so that both effects seem to push in the same direction [jl^, [l9| ]. 

4. Final remarks 

Strong magnetic fields can modify the nature of the chiral and the deconfining transitions, 
opening new possibilities in the study of the phase diagram of QCD, introducing a new control 
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parameter, besides temperature and baryon chemical potential, in the study of the thermodynamics 
of strong interactions. These high magnetic fields are also essential in the context of high-energy 
heavy ion collisions in generating a measurable charge asymmetry to determine the presence of 
CP-odd domains created by sphaleron transitions. Several relevant questions that are still not fully 
covered can be raised. How strong can one make these magnetic fields in current and future exper- 
iments? How long lived could they be, and how uniform, for different values of y/s? Only the first 
theoretical estimates [§, 42] and preliminary lattice simulations [13] were performed, providing 
encouraging results. 

On the one hand, due to the presence of a strong magnetic background, non-central heavy 
ion collisions might show features of a first-order transition when contrasted to central collisions. 
On the other hand, finite-size effects are sizable for non-central collisions [p^], so that an accurate 
centrality dependence study seems to be necessary. This demands a thin binning of centrality and 
control of finite-size effects, a very difficult (but necessary) task for experimentalists in the process 
of data analysis, especially in the scheduled Beam Energy Scan program at RHIC-BNL. 

Finally, one still needs to perform dynamical investigations within the chiral magnetic effect 
scenario to determine the relevant time scales and verify whether effects from the CP-odd domains 
survive long enough to produce measurable signatures. 
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